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Abstract 

In this paper we derive an effective equation for derivative pricing which 
accounts for the presence of virtual arbitrage opportunities and their ehmination 
by the market. We model the arbitrage return by a stochastic process and find an 
equation for the average derivative price. This is an integro-differential equation 
which, in the absence of the virtual arbitrage or for an infinitely fast market 
reaction, reduces to the Black-Scholes equation. Explicit formulas are obtained 
for European call and put vanilla options. 



1 Introduction 

The Black-Scholes (BS) analysis |I|] of derivative pricing is one of the most beautiful 
results in financial economics. The BS formulas are easy to understand and to handle 
which leads to wide use of the formulas by traders. There are several assumptions 
in the basis of BS analysis such as the quasi-Brownian character of the underlying 
price process, constant volatility and, the most importantly for our goal, the absence 
of arbitrage. Though first two assumptions can be relaxed and a lot of attempts to 
improve BS analysis using stochastic volatility models and price models beyond the 
quasi-Brownian class have been made, the third assumption, namely the no-arbitrage 
condition has not been tackled yet. Indeed, almost all models imply the no-arbitrage 
constraint, employ the martingales and hedging portfolios to calculate derivative prices. 
More generally, the no-arbitrage constraint is at the heart of contemporary financial 
calculus 0, ^ and leads to many fruitful results not only in derivative pricing but also 
in portfolio theory and asset pricing [||, ||. 

Accepting that the local (virtual) arbitrage is short lived and hence irrelevant for 
long term pricing, one is forced to consider the arbitrage if the study of the short term 
behavior is one's objective. It is well-known that local arbitrage opportunities are 
always present in the market |^ and have finite life-time (for example, for extremely 
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liquid markets such as futures on S&P the characteristic arbitrage time is of the order 
of 5 minutes |^ ; for less liquid markets such as the bonds market the time can be much 
longer). This poses the question of how to account the local arbitrage opportunities 
and improve the pricing formulas. 

One of the possible approaches to the problem has been suggested in Ref 0. It was 
argued there that the violation of the no-arbitrage constraint and the non-Brownian 
character of underlying price walks can be accounted for in the same framework. Being 
a consistent theory, the approach is however extremely complicated and does not allow 
one yet to get simple analytical results which would be easy understandable and han- 
dleable. That is why in this paper we develop simplified tractable analytical version to 
account for virtual arbitrage opportunities which does not use complicated techniques 
and results in simple enough final equations. 

The paper is organized as follows. In next section we formulate the model. To this 
end the BS equation is rederived and it is shown how the local arbitrage opportunities 
can be introduced in the model. In section 3 we find an equation for the average 
derivative price. This equation generalizes the BS equation and converges to it in 
the limit of absence of the arbitrage opportunities or infinitely fast market reaction. In 
conclusion we discuss drawbacks of the model and possible ways to improve it. Explicit 
analytical solutions for European vanilla options are contained in Appendix. 

2 Formulation of the model 

We start this section with a quick standard derivation of the BS equation. This both 
makes the paper self-contained, introduces useful notation, and allows us to emphasize 
the important notion which is extensively used later. 

First, let us denote V(t, S) as the price of a derivative at time t condition to the 
underlying asset price equal to S. We assume that the underlying asset price follows 
the geometrical Brownian motion, i.e. 

— = fidt + adW , 

with some average return /i and the volatility a. They can be kept constant or be 
arbitrary functions of S and t. The symbol dW stands for the standard Wiener process. 
To price the derivative one forms a portfolio which consists of the derivative and A 
units of the underlying asset so that the price of the portfolio is equal to 11: 

The change in the portfolio price during a time step dt can be written as 

... ..r ..n fdV a^S^d^V\ ^ (dV ,\ ^ 
dn^dr-Ad^^ - + ^^ dt+ --A d^. 



from of Ito's lemma. We can now chose the number of the underlying asset units 

dV_ 

as 



A to be equal to -^ to cancel the second term on the right hand side of the last 



equation. Since, after cancellation, there are no risky contributions (i.e. there is no 
term proportional to dS) the portfolio is risk-free and hence, in the absence of the 
arbitrage, its price will grow with the risk-free interest rate r: 

dn = rUdt , (1) 

or, in other words, the price of the derivative V{t, S) shall obey the Black-Scholes 
equation: 

In what follows we use this equation in the following operator form: 

C^sV = , ^BS=g-^ + ^g^ + rS--r. (3) 

To formulate the model we return back to Eqn(|l]). Let us imagine that at some moment 
of time T < t a. fluctuation of the return (an arbitrage opportunity) appeared in the 
market. It happened when the price of the underlying stock was S' = S{t). We then 
denote this instantaneous arbitrage return as z/(r, S"). Arbitragers would react to this 
circumstance and act in such a way that the arbitrage gradually disappears and the 
market returns to its equilibrium state, i.e. the absence of the arbitrage. For small 
enough fluctuations it is natural to assume that the arbitrage return TZ (in absence of 
other fluctuations) evolves according to the following equation: 

^ = -A7^ , 7^(r) = u{t, S') (4) 

at 

with some parameter A which is characteristic for the market. This parameter can 
be either estimated from a microscopic theory like |^ or can be found from the mar- 
ket using an analogue of the fluctuation-dissipation theorem 0. In the last case the 
parameter A can be estimated from the market data as 



TTlog 



(t - t') 



C^BsV ,,, CbsV 



W7^^^W(^') 
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and may well be a function of time and the price of the underlying asset. In what follows 
we however consider A as a constant to get simple analytical formulas for derivative 
prices. The generalization to the case of time-dependent parameters is straightforward. 
The solution of Eqn(^ gives us TZ(t, S) = z/(r, 5')e~'^^*~'^'' which, after summing over 
all possible fluctuations with the corresponding frequencies, leads us to the following 
expression for the arbitrage return at time t: 

ft /"OO 

7^(t,5)=/dr/ dS' e-^^'-^^Pit,S\T,S')i^iT,S') , t<T (5) 

JO Jo 

where T is the expiration date for the derivative contract started at time t = and the 
function P{t, S\t, S') is the conditional probability for the underlying price. To specify 



the stochastic process z/(t, 5*) we assume that the fluctuations at different times and 
underlying prices are independent and form the white noise with a variance T? -fit): 

{u{t,S)) = , {u{t,S)u{t',S')) = J:^-e{T-t)f{t)6{t-t')6{S -S') . (6) 

The function f{t) is introduced here to smooth out the transition to the zero virtual 
arbitrage at the expiration date. The quantity S^ ■ f{t) can be estimated from the 
market data as 

and has to vanish as time tends to the expiration date. 

Since we introduced the stochastic arbitrage return TZit, S), Eqn([^) has to be sub- 
stituted with the following equation: 

dn = [r + 7^(t, S)]Udt , (7) 

which can be rewritten as 

CBsV = nt,S)(v-S^] , (8) 

using the operator C^s- Eqns(|^),(^ and (|^) complete the formulation of the model. 

It is worth noting that the model reduces to the pure BS analysis in the case of 
infinitely fast market reaction, i.e. A —>■ oo. It also returns to the BS model when there 
are no arbitrage opportunities at all, i.e. when S = 0. 

In the presence of the random arbitrage fluctuations TZ{t, S), the only objects which 
can be calculated are the average value and other higher moments of the derivative 
price. In this paper we examine the average price and derive the pricing equation for 
it in the next section. 



3 Effective equation for derivative price 

We start this section with the note that the probability distribution of TZ(t, S) is Gaus- 
sian. Moreover, it can be shown that the probability (up to a normalization constant) 
of the trajectory Tl{-, ■) has the form: 



P[7^(•,■)]~exp 



dtdt'd^d^' 7^(t, s)K-\t, s\t', s')n{t', s') 



(9) 



2S2 7o 
where the kernel of the operator K is defined as: 

/■oo 

K{t, S\t\ S') = e{T - t)e{T - t') / drds f{T)e{t - T)e{t' - r)e-^(*+*'^'^) 

JO 

xP(t,5|r,s)P(t',5'|r,s) . (10) 

It is easy to see that the kernel is of order 1/A and vanishes as A ^ cxd. Eqn(P), in 
particular, results in the equality for the correlation function: 

(7^(t, s)n{t', s')) = T? ■ K{t, s\t', s') , (ii) 



which we use below. 

Now let us return to the dynamical equation for the derivative price 

>CBs^(t, S) = n{t, S) [Vit, S) - SdsV{t, S)] (12) 

and note that, since S^/A plays a role of small parameter in the problem, the noise TZ 
can be considered as weak and we can find a formal iterative 7?.-dependent solution of 
the last equation. In the lowest non-trivial order we have the equation: 

/:bs^ = n[V - SdsV] = 7^[l - Sds]CBs''Tl[V - SdsV] , (13) 

which after averaging (using (0)) over all possible realizations of the fluctuations TZ 
give us an equation for the average derivative price V = {V)ti up to and including 
terms proportional to T? /\: 

oo 

CBsVit, S) = S^ fdt'dS' [(1 - 595)/:bs] (t^ S\t', S')-K{t, S\t', S') [(1 - Sds)v] {t', S') 



(14) 
together with the payoff condition: 

y(T,5) = Vpayoff(5). 



Equation (|T^) is the central result of the paper. This is an integro-differential 
equation which in the limit A — >• cx) or S — i> reduces to the Black-Scholes equation 
and accounts for local arbitrage opportunities and the corresponding market reaction 
to them. Equations of this type are very familiar in physics where they are called 
one-loop effective Dyson-type equations. 

To conclude the section, it is interesting to note that due to properties of the 
integrand on the right hand side of Eqn(|T^, the integration is effectively limited to 
the interval from time t to the expiration date. It means that any mispricings which 
happened in the past do not influence the derivative price, as one would expect, and 
the only relevant contribution comes from future mispricings. The explicit solutions 
for European vanilla options are derived in the Appendix. 

4 Conclusion 

In conclusion we want to discuss some obvious drawbacks of the model and ways to 
improve it. 

First of all, all critical comments of BS analysis can be forwarded to this model, 
except for the no- arbitrage constraint. Indeed if there are transaction costs or if the 
price process for the underlying asset is not quasi-Brownian motion, it is impossible to 
create a risk-less portfolio and, hence, to derive the model. This is not a new problem 
and many efforts to overcome this difficulty have been undertaken. It is possible to 
demonstrate that the virtual arbitrage model can be improved in the same manner by 
these methods as they succeed for no-arbitrage BS analysis. This, in principle, allows 
one to include the transaction costs, the "fat" tails for the underlying assets probability 



distribution function and the stochastic volatihty in the present model by redefinition 
of the operator £bs- 

The second point concerns the market reaction to the arbitrage opportunity, or, 
quahtatively the form of TZ(t, S) in Eqn(|^). It may be argued that the market reaction 
is not exponential as assumed in Eqn(||), but has another functional dependence. This 
dependence can be found from statistical analysis of the stock and derivative prices and 
then included in the equation for TZ{t, S) (in particular, the functional dependence can 
change with time, for example A in (H) can be a function of r). It certainly complicates 
the model but leaves the general framework intact. 

Another point to consider is the absence of correlations between virtual arbitrage 
opportunities which we assumed in the text, i.e. the white noise character of the 
process z/(t, S). It is clear that some correlations can be easily included in the model 
by substituting the relations in Eqn(^ by the equations: 

{u{t, S)) = , {u{t, S)u{t', S')) = T? ■ Fit, S\t', S') 

with some correlation function F{t,S\t' ,S'). Such generalization, though making the 
analytical study almost impossible, allows one to proceed with numerical analysis for 
the model. 

Finally, the model contains new parameters such as S and A. Though the param- 
eters can be estimated from the statistical data as we mentioned above, they still can 
be considered as a fitting parameters. If we follow this line, to define their values some 
kind of test of the final formulas on real market data should be carried out. The least 
curve implied volatility can be used as an example of such tests. 
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Appendix: Solution for European vanilla options 

Before going further we would like to point out that, since Eqn(|lJ) is linear, there 
exists the call-put parity theorem for the average prices C(t, S", T, E) and P(t, S", T, E) 
of European call and put with the same strike price E and the same expiration date 
T: 

C(t, S, T, E) - Pit, S, T, E) = Fit, S, T, E) , 

where F(t, S, T, E) is the price of the corresponding forward agreement. It means that 
pricing of call and put also gives one a price for the forward agreement. 

In this appendix we derive explicit formulas for European call and put options for 
the case of geometrical random walk for the price of the underlying asset. In this 



case it is convenient to use the variable y = \og{S/E) instead of S. The probabihty 
distribution function then takes the form: 



a^2Tr{t - t') 
and the operator £bs can be rewritten as: 



exp 



l\\2 



{y-y'-^i{t- t')} 
2(r^{t-t') 



CBS = dt + Y9l+(r-Y^d,-r. 



^y 



The corresponding kernel K{t, y\t', y') is equal to 

/•OO 

K{tMt\y) = e{T-t)e{T-t')j^ dr /(r^t - r^t' - r)e-^(*+* -2-) 

{y-y'-^j,{t-t')Y' 



aJ27r{t + t' -2t) 



exp 



2(r^{t + t' -2r) 



(15) 



In the limit of sufficiently fast market relaxation (A ^ cr^, /i) and the particular 
choice of the smoothing function f{t) 



the kernel takes the form 

K{t, y\t', y') = 5{y - y')9{T - t)9{T - t')K{t, t') 
where the following notation has been introduced 



(16) 



(17) 



K{t,t') = ^^(t - t') |e-^(*^*') 
+ ^^(t' - t) |e-^(*'-*) 



1 
1- 



Ae-27(T-t') 

A + 7 
Ae-27(T-t)- 

A + 7 



-X{t+t') 



_ Q-Ht+t') 



1- 
1- 



Ae^^ 

A + 7 
Ae^^ 

A + 7 



If A(t + t') ^ 1, that is t, t' are sufficiently far from t = 0, then the terms proportional 
to e"'^*^*"'"* ) in (|l^) can be neglected and we obtain for K(t,t'): 



K{t,t') = ^e{t-t')e-'^'-''^ 



1 



Ae-^^(^-*') 
A + 7 



+ ^^e{t' - t)e-^(*'-) 



Ae-27{T-t) 



A + 7 



A kernel of the integral operator £gg in this case can be written as 



^Bsit^y\t'^y' 



e{t' - 1) 



oJ2Ti{t' - 1) 



:exp 



{y'-y-a{t'-t)f 
2a2(t' - 1) 



r{t' - 1) 



(19) 
(20) 



with the parameter a given by the relation 



a 



a = r — 



f211 



Now we will solve Eqn ([l^ keeping in mind that T? /\ is a small parameter in the 
problem. At the first non-trivial order we have: 

V = V, + ^-V^ (22) 

where the zero-order solution Vo(t, y) is the BS solution for the derivative: 

C^sVoit, y) = , \/o(T, y) = l^payoff (y) 
and Vi satisfies the equation 

oo oo 

CBsHt.y) = Jdt' J dy' [(1 - dy)C^l\ {t,y\t',y')K{t,y\t',y') [(1 - dy)Vo] {t',y') (23) 

-oo 



which can be explicitly rewritten as: 

CBsVi{t,y) = F{t,y) , Vi{T,y) = 
with the right hand side given by the equation: 



(24) 



Fit,y) = -^(^ + ^]0{T-t) 



4a 



cr" 



A 



A + 7 



-27(r-i) 



X J^ dt' ^—^==^ [(1 - dy)Vo] {f, y) . 



271(1' - 1) 



(25) 



Following these general formulas, the prices of the European call and put options 
have the form: 

C(t, y) = Coit, y) + \Ciit, y) , P{t, y) = P^{t, y) + ^Pi(t, y) . 

The zeroth order solutions are just BS solutions for the options: 

Co{t,y) = E[e^iV[di(y,T-t)]-e-''(^-*)iV[d2(l/,T-t)]] (26) 

Poit^y) = E\e-^'^^-'^N[-d2{y,T -t)]-eyN[~d,iy,T -t)]\ (27) 



with the notation 



di{y,r) 



y + {a + a^)T 



(JxT 



d2{y,T) 



y + ar 



and 



N[x) 



dy e 



-J/V2 



The correction terms Ci{t,y), Pi{t,y) can be obtained after some straightforward al- 
gebra: 



Ci{t,y) 



Ee-^'^-t) ( 1r 



T-t 



4a 



o' 



dr {T-t-T + 







e-27(r-t) _ g-27T ■ 

27(1 + 7/A) 



X ^= N[d2{y, T-t-T)] 



Pi{t,y) 



V27rr 

Ee-^'(^-^) 
4^ 



(28) 



2r 



a" 



T~t 



dr <^T-t-r + 



g-27(T-t) _ g-27r 

27(1 + 7/A) 



-(A+aV2fT^)^ 



X- 



V2 



TTT 



-N[-d2iy,T-t-T)] 



(29) 



When the expiration occurs at sufficiently large time that we can neglect its influ- 
ence, the smoothing function f(t) can be substituted by 1, i.e. the limit 7 — > 00 can 
be taken. In this case expressions (^,29) have an even simpler form: 



Ci{t,y) 



Pi{t,y) 



^e-^(^-*) / ^ 2rN^ 



4a 



Ee-'-(^-*) 
4^ 



a 



2 / 



T-t 

f dr (T-t-T) 



-iX+a'^/2a'^)T 



T-t 



V2nT 



N[d2iy,T-t-T)] (30) 



1 + - dr{T-t-r) ^^N[-d2{y,T-t-T)]. 



(31) 
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